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摘 要 : 时 间 频 率 数据 资源 是 国家 重要 的 信息 资源 ， 时 间 频 率 数 据 平台 设计 依托 国家 授 
时 中 心 运行 产生 的 时 频数 据 、 国 际 权 度 局 (Bureau International des Poids et Mesures , 
BIPM) 、 国 际 地 球 自 转 服务 机 构 (International Earth Rotation and Reference Systems 
Service , IERS) 公布 的 时 间 频 率 相关 数 据 等 实现 。 该 数据 平台 可 为 守 时 技术 研究 提供 丰 
富 的 数据 资源 支撑 ,首先 介绍 时 间 频 率 数据 平台 的 建设 情况 , 随后 利用 数据 平台 汇 交 的 本 地 
守 时 系统 比 对 、 远 距离 链 路 比 对 数据 开展 守 时 技术 研究 ， 包括 原子 钟 状态 评估 ,主要 分 析 原 
子 钟 异 常 跳 变 情况 ， 研 究 跳 变数 据 处 理 方 法 ， 针 对 不 同类 型 守 时 钟 ， 分 析 对 比 波 动情 况 、 稳 
定 度 等 原子 钟 性 能 指标 ， 此 外 基于 ARIMA 模型 开展 钟 差 预 报 研究 。 以 上 研究 结果 可 为 守 时 
系统 连续 可 靠 运行 提供 重要 参考 。 
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随 着 我 国 科技 不 断 进步 ， 时 间 已 经 成 为 非常 重要 的 研究 课题 趾 。 高 精度 时 间 对 于 各 行 各 
业 的 影响 巨大 ， 包 括 卫星 导航 、 计 量 等 科技 研究 方面 ， 交 通 运输 、 电 力 配 送 、 社 会 生产 等 生 
活 多 方面 "站 。 对 于 时 间 的 研究 ， 在 远 距 离 比 对 方面 ， 基 于 北斗 系统 的 高 精度 时 间 比 对 技术 
研究 前 景 广阔 ， 在 本 地 守 时 系统 方面 ， 开 展 时 间 保 持 技 术 的 研究 、 原 子 钟 物理 性 能 研究 、 守 
时 系统 驾驭 控制 相关 研究 等 意义 非凡 ， 这 一 切 研究 的 基础 依托 于 时 间 频 率 数据 资源 的 汇 交 。 

对 日 益 增 长 的 时 间 频 率 数据 资源 服务 需求 ， 时 间 频 率 数据 平台 应 运 而 生 广 1。 该 数据 
平台 的 建立 在 时 间 频 率 数据 规范 化 存储 、 数 据 共享 、 数 据 网 络 形成 、 用 户 基础 服务 、 科 研 驱 
动 等 方面 具有 重要 作用 。 数 据 平 台 采 用 多 层 应 用 结构 组 件 化 开发 模式 ， 数 据 分 类 定 级 ， 整 个 
系统 稳定 性 强 ， 能 够 满足 一 定 规模 的 应 用 需求 ， 具 有 高 效 、 安 全 、 可 扩展 等 特点 。 目 前 建成 
的 数据 服务 平台 汇集 了 本 地 守 时 系统 、 远 程 链 路 比 对 、 国 际 服 务 机 构 发 布 相关 数据 ,涵盖 了 
原子 钟 权重 、 原 子 钟 速率 、 原 子 钟 频率 漂移 、CGGTTS、RINEX、 公 报 数据 等 ， 在 数据 结构 
层面 ， 较 为 全 面 地 对 数据 进行 分 级 分 类 ， 基 本 满足 用 户 一 站 式 应 用 需求 。 

用 户 可 以 依据 时 间 频 率 数据 平台 对 时 频数 据 资 源 进行 深度 挖 掘 ， 开 展 各 项 基础 研究 ， 
发 挥 其 应 有 价值 。 本 文 主 要 从 原子 钟 性 能 分 析 、 原 子 钟 异常 跳 变 、 钟 差 预报 展开 ， 这 三 个 方 
面 是 守 时 系统 稳定 运转 的 重要 基础 研究 参考 信息 , 对 于 保障 守 时 系统 连续 可 靠 运转 具有 重要 
意义 。 当 然 ， 时 频数 据 的 应 用 不 拘 于 此 ， 未 来 也 将 汇集 更 多 相关 数据 资源 ， 开 展 更 多 层次 、 
多 方位 研究 。 
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时 间 频 率 数据 对 天 文学 、 地 球 动力 学 、 物 理学 研究 ， 以 及 信息 传递 、 电 力 输 配 、 深 空 探 


测 、 空 间 旅行 、 导 航 定位 、 武 器 实验 、 地 震 监 测 、 计 量 测试 等 方面 有 重要 的 应 用 价值 ， 因 此 
建设 专属 时 间 频 率 数据 平台 意义 深远 中 。 
时 间 频率 数据 平台 依据 时 间 频 率 数据 结构 特征 进行 规范 化 存储 , 整个 数据 平台 依据 多 层 


应 用 结构 设计 ， 数据 平台 


安全 性 。 


E 最 核心 也 是 最 底层 部 分 ,不 会 直接 暴露 ， 可 以 有 效 提高 数据 库 的 


时 间 频 率 数 据 平 台数 据 模块 划分 按时 间 频 率 数 据 分 类 定 级 标准 完成 。 分 类 定 级 的 标准 包 
括 数据 的 使 用 方法 、 数 据 开 放 计划 、 数 据 本 身 的 特性 、 数 据 的 更 新 频率 、 数 据 的 体 量 大 小 等 
进行 划 定 。 目 前 平台 的 数据 资源 包括 原子 钟 比 对 数据 、 原 子 钟 状态 参数 数据 、 国 际 权 度 局 发 


布 的 数据 、 公 报 数据 等 。 


系统 平台 的 整体 扩展 性 较 强 ， 随 着 原子 钟 的 增加 、 退 出 、 守 时 系统 模块 的 更 新 , 访问 的 
需求 量 增加 , 可 按照 新 的 需求 通过 在 不 同系 统 层面 调度 升级 组 件 或 扩展 新 加 入 的 组 件 来 升级 
系统 ， 以 适应 新 的 需求 及 不 断 变化 的 数据 库 建设 要 求 。 


建成 的 数据 平台 具有 安全 性 、 


完整 性 、 可 靠 性 、 连 续 性 和 权威 性 ， 可 以 为 时 频 科学 研究 


与 发 展 提供 重要 的 基础 技术 支撑 , 同时 也 为 各 项 时 频 技术 和 应 用 研究 提供 基本 保障 ,图 1 为 


时 间 频 率 数据 平台 主页 。 


Ji 时 间 频 率 数据 平台 


Fig.1 Time frequency data platform 


2 基于 数据 平台 相关 数据 分 析 研 究 


2.1 原子 钟 异 常 跳 变 及 处 理 方 法 


基于 时 间 频 率 数 据 平 台 本 地 守 时 系统 测 


子 钟 异常 跳 变 及 相关 处 型 


方法 研究 。 


数据 、 原 子 钟 指标 实时 测量 数据 , 可 以 进行 原 


la 


原子 钟 跳 变 分 为 相位 跳 变 和 频率 跳 变 0， 其 中 相位 跳 变 较为 明显 ， 对 于 整个 守 时 系统 的 
扰动 相对 强烈 ~“， 故 本 文 主要 讨论 相位 跳 变 。 相 位 跳 变 的 原因 有 很 多 ， 主 要 包括 原子 钟 本 


身 物理 异常 、 所 处 环境 温度 、 湿 度 变化 、 电 磁场 干扰 和 电源 中 断 等 。 异常 跳 变 对 于 守 时 系统 


的 稳定 运行 有 一 定 影响 ， 例 如 原子 钟 跳 变 会 使 得 连续 稳定 运行 的 守 时 系统 产生 波动 ， 此 外 ， 
如 果 该 原子 钟 未 被 限制 权重 ， 则 导致 系统 计算 的 要 驭 量 异 常 。 
异常 跳 变 检测 是 基于 门限 阔 值 , 超出 门限 即 判定 为 突变 , 检测 门限 值 往往 根据 经 验 值 设 
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定 , 检 出 异常 值 后 ， 处 理 异常 跳 变 的 方法 有 最 小 二 乘 拟 合法 ， 即 计算 最 小 偏差 平方 和 得 到 拟 


合 多 项 式 ， 原 理 为 假定 拟 合 多 项 式 


y=agptaxt+L tan” ， (1) 


则 拟 合 偏差 平方 和 表示 为 


e= db, —(a,+a,x+L +a,x")P ; (2) 


i=l 


最 小 值 即 偏 导数 为 0 求 出 的 极 值 ， 针 对 上 式 右边 求偶 


m 
1 x L x |e, 


令 
1 为 工 x a, 
xa L X A= a, 
MMO M M 
lx Lx Ay 
则 拟 合 值 为 


A=(X'X)(X'Y) 


1 x, L x ol 
MMO M|M 


1 x, L x an 


导数 并 化 简 可 得 


yı 
X2 (3) 
Yn 
vy 
Y2 
， (4) 
Vn 
(5) 


从 时 频数 据 平 台 随 机 选择 中 国 科 学 院 国 家 授时 中 心 一 台 钨 原子 钟 Cs3436， 进 行 扣 常 数 
转 纳 秒 ， 假 定 在 儒 略 日 59370.253 处 发 生 相 位 跳 变 , 如 图 2 。 运 用 最 小 二 乘 拟 合法 对 跳 变 进 
行 修正 ， 结 果 如 图 3。 由 图 3 可 知 ， 该 方法 对 于 异常 值 跳 变 修正 情况 良好 ， 证 明 算法 有 效 。 


x107 
> 


clock difference/ns 
is: 
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图 2 Cs3436 钟 差 波 动情 况 


59390 59395 


Fig.2 The clock difference fluctuation of Cs3436 
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图 3 原子 钟 异常 跳 变 处 理 结果 


Fig.3 Processing result of abnormal jump of atomic clock 


2.2 守 时 型 原子 钟 性 能 分 析 


利用 时 间 频 率 数 据 平台 公报 数据 、 原 子 钟 相位 数据 、 速 率 数据 、 远 距离 链 路 比 对 数据 等 
开展 守 时 型 原子 钟 的 性 能 分 析 "" 。 稳 定 度 是 评估 原子 钟 的 重要 指标 ， 艾 伦 (Allan) 方差 
是 时 间 领 域 常用 的 分 析 时 间 频 率 稳 定 度 的 测量 方法 , 具有 良好 的 置信 度 。 艾 伦 方差 计算 公式 
为 


1 N-2 

2 2 

CT = 一 一 一 一 一 》 X — 2X., +X, 6 
i M HD A mt) < 


>» 


Hop, 为 第 ' 个 相位 数据 间隔 ，ML 为 连续 测量 的 次 数 ，T 为 采样 间隔 。 从 时 间 频 率 
数据 平台 随机 选择 中 国 科学 院 国 家 授时 中 心 3 台 钨 钟 和 3 Gaver 2022 年 1 月 份 比 对 数据 ， 
相关 波动 情况 如 图 4， 稳定 度 计算 结果 如 图 5。 
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4 氧 原 子 钟 和 钨 原子 钟 波动 情况 对 比 
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Fig.4 Comparison of fluctuations of hydrogen maser and cesium atomic clock 
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图 5 氧 原子 钟 和 钨 原子 钟 稳定 度 对 比 


Fig.5 Comparison of the stability of hydrogen maser and cesium atomic clock 
图 4 数据 处 理 的 过 程 为 选取 一 个 月 原子 钟 钟 差 数 据 进 行 速率 及 频率 漂移 预报 , 随后 进行 
频率 修正 ,由 图 4 可 以 看 出 , 钨 原子 钟 钟 差 波 动 变化 较为 明显 , 氧 原 子 钟 钟 差 波动 变化 较 小 ， 
但 是 长 期 具有 趋势 项 。 图 5 为 利用 艾 伦 方差 计算 6 台 钟 稳定 度 对 比 , 可 以 看 出 氨 原 子 钟 短期 
稳定 度 高 于 铭 原 子 钟 ， 但 由 于 氢 原 子 钟 普遍 频率 漂移 较 大 ， 长 期 稳定 度 变 差 。 


2. 3 基于 ARIMA 模型 的 钟 差 预报 


原子 钟 的 钟 差 预报 是 守 时 系统 动态 控制 参考 量 计 算 的 基础 "“ ， 同 时 也 有 原子 钟 异 常 检 
测 的 功能 ， 当 预报 钟 差 与 实际 钟 差 相差 大 于 一 定 范围 , 则 判定 为 原子 钟 异 常 ， 因 此 钟 差 的 预 
报 具有 重要 意义 上 。 由 于 原子 钟 钟 差 数据 往往 是 非 平 稳 序列 ， 而 ARIMA 模型 将 非 平 稳 数 
据 经 过 差分 得 到 平稳 数据 ， 建 模 简单 ， 运 算 过 程 只 依赖 数据 本 身 ， 故 本 文采 用 ARIMA 模型 
进行 钟 差 预报 ，ARIMA 模型 由 AR 模型 及 MA 模型 两 部 分 构成 ，AR 模型 即 自 回归 模型 ， 
可 以 表示 为 


X=hx +hx ,+L +h 


pXp-2 TU, (7) 


? 


其 中 , P 为 自 回 归 模 型 的 阶 数 ; bh, AA 


a 


归 系 数 。 通常 认为 AR 模型 中 , |e g 阶 移动 平均 。 
MA 模型 表示 为 


u=E,+PE,+BE,E +QE (8) 


其 中 ，8 NAMB. ARIMA 模型 一 般 表达 式 为 


X=Ax thekb +h yt @&_,,~ 6 t+ Pio (9) 


ARIMA 模型 的 核心 是 通过 历史 数据 预测 当前 数据 ， 对 于 非 平 稳 数 据 ， 进 行 差 分 获得 平 
稳 数据 作为 模型 输入 。 模 型 的 自 回 归 阶 数 P ,滑动 平均 阶 数 4 以 及 差分 阶 数 d 需要 确定 ， 其 
中 差分 阶 数 d 的 确定 是 数据 经 过 差分 成 为 平稳 序列 的 差分 次 数 ， 自 回归 阶 数 了 PD 、 滑 动 平均 
阶 数 4 可 以 通过 求 取 自 相关 函数 和 偏 自 相关 函数 获得 ， 也 可 以 基于 AIC 准则 、BIC 准则 获 


得 


= 


基于 时 频数 据 平台 随机 选择 中 国 科 学 院 国 家 授时 中 心 一 台 氢 钟 VCH067 一 个 月 的 数据 
作为 研究 基础 ， 运 用 ARIMA 模型 进行 预报 ， 钟 差 预报 的 结果 如 图 6， 钟 差 预报 与 实际 钟 差 
RAIA 7， 钟 差 预 报 的 RMS 如 图 8。 
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图 6 钟 差 预报 结果 
Fig.6 Clock difference forecast result 
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图 7 钟 差 预 报 误差 
Fig.7 The prediction error of clock difference forecast 
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图 8 钟 差 预报 的 RMS 


Fig.8 The prediction RMSE of clock difference forecast 


由 图 6 和 图 7 可 以 看 出 ， 基 于 ARIMA 模型 钟 差 预报 结果 与 实际 钟 差 数据 吻合 度 恨 好 ， 


趋势 一 致 。 预 报 钟 差 与 实际 钟 差 的 残 差 保持 在 土 0.2ns 以 内 。 图 8 是 以 天 为 滑动 区 间 长 度 计 
算得 到 的 预报 钟 差 的 均 方 根 ， 均 方 根 基本 保持 在 0.1ns 以 内 ， 说 明基 于 ARMA 模型 进行 的 
钟 差 预报 效果 良好 。 该 结果 可 为 守 时 系统 的 钟 差 预 报 、 卫 星 钟 差 预 报 提供 一 定 参考 。 


3 结语 


本 文 首先 介绍 时 频数 据 平台 建设 背景 及 具体 部 署 情况 , 随后 针对 数据 平台 本 地 原子 时 比 


对 数据 、 国际 公报 等 数据 开展 基础 研究 ， 相 关 研 究 包 括 原子 钟 的 中 变 分 析 及 处 理 方法 、 守 时 
型 原子 钟 波动 情况 、 稳 定 度 性 能 分 析 ， 以 及 基于 ARIMA 模型 的 钟 差 预 报 研 究 。 三 种 基础 研 
究 可 为 守 时 系统 稳定 可 靠 、 连 续 运 转 提 供 参 考 。 未 来 随 着 时 间 频 率 数 据 平台 建设 的 不 断 完善 、 
数据 集 种 类 的 不 断 丰 富 、 优 化 ， 其 应 用 空间 必 将 更 加 广阔 。 
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Abstract: Time-frequency data resources are important national information resources. The 
design of the time-frequency data platform relies on the time-frequency data generated by the 
daily operation of the National Time Service Center, the bulletin issued by Bureau International 
Des Poids et Mesures (BIPM) and International Earth Rotation and Reference Systems 
Service(IERS).The data platform can provide abundant data resources for the research of 
time-keeping technology. This paper first introduces the construction of the time-frequency data 
platform, and then uses the local time-keeping system comparison and long-distance link 
comparison data collected by the data platform to carry out time-keeping technology research, 
including the evaluation of the state of the atomic clock, and the analysis of the abnormal jump of 
the atomic clock. For different types of clock keeping, analyze and compare the performance 
indicators of atomic clocks such as fluctuation and stability, and conduct clock error prediction 
research based on the ARIMA model. The research contents of those different aspects can provide 
important reference for the continuous and reliable operation of the time-keeping system. 

Key words: time frequency; data platform; time-keeping technology; atomic clock; time-keeping 


system 


